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Abstract The spin up of helium II is studied by calculating the spin-down 
recovery of a superfluid-filled container after an impulsive acceleration and 
' comparing with experiments. The calculation takes advantage of a recently 

published analytic solution for the spin up of a Hall-Vinen-Bekharevich- 
Khalatnikov superfluid that treats the back-reaction torque exerted by the 
viscous component self-consistently in arbitrary geometry for the first time. 
Excellent agreement at the 0.5% level is obtained for experiments at T = 
1.57 K, after correcting for the non-uniform rotation in the initial state, con- 
firming that vortex tension and pinning (which are omitted from the theory) 
play a minimal role under certain conditions (small Rossby number, smooth 
walls) . The dependence of the spin-down time on temperature and the mass 
fraction of the viscous component are also investigated. Closer to the lambda 
point, the predicted onset of turbulence invalidates the linear Ekman theory. 

Keywords Rotation • Superfluid • Helium II 



1 Introduction 

o 



The spin up of helium II in closed vessels has been studied experimentally 
[TTl[30lfT^[15l[20l[23l[24l[25] and theoretically [4l[T][T6|[T9] . but good agreement 
has proved elusive. The problem poses a significant challenge over classical 
Navier-Stokes spin up because of the additional complications arising from 
the two-component hydrodynamics (coupled by mutual friction) in a super- 
fluid and quantum mechanical vortex pinning and nucleation. Early studies 
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of helium II spin up measured the torque required to accelerate a fluid- 
filled container, while monitoring the fluid response with a Rayleigh disk 
device suspended in the flow [11,30,12,15 . Asymmetric acceleration and re- 
tardation and strong hysteresis were observed, effects which are enhanced in 
rough-walled containers, suggesting that vortex pinning and nucleation play 
a crucial role. 

The first comprehensive experimental study of the spin up of helium II 
was conducted by Tsakadze and Tsakadze 23,24,25 , who monitored the 
hydrodynamic relaxation of superfluid-filled containers after an impulsive 
increase in angular velocity. Spherical and cylindrical containers made of 
plcxiglas or duralumin were used; an additional coating of plexiglas powder 
was applied to study strong vortex pinning. For smooth- walled containers, a 
clear decrease in the spin-up time with increasing temperature is observed. 
The spin-up time is also a function of initial and final angular velocity, con- 
tainer size, and density fraction of the viscous component. In rough-walled 
containers, a significantly shorter relaxation time is observed, which does not 
depend on the temperature or velocity jump. The above investigations have 
been repeated in spherical Couette geometry [25] . 

No complete theoretical explanation of the Tsakadze results has been 
published to date. The lack of a theory is felt most keenly in the field of neu- 
tron star astrophysics, where the Tsakadze experiments serve as a terrestrial 
facsimile of the glitch phenomenon. The origin of glitches is a fundamental 
40-year-old puzzle [27] which relates to the hypothesis of nuclear superfluid- 
ity [27] . Pioneering attempts to theoretically understand the nature of spin 
up in helium II focused on the interaction between quantized vortices in 
the superfluid and the boundaries [U[T] ■ A rigorous analytical treatment was 
provided by Reisenegger |19j . who solved the vortex- fluid equations of mo- 
tion of Chandler and Baym [5] between two infinite, rotating parallel plates 
accelerating at a slow and constant rate. In smooth-walled containers, two 
spin-up mechanisms operate; classic Ekman pumping in the viscous compo- 
nent near the superfluid transition temperature, and Ekman-like circulation 
in the inviscid component near absolute zero. In the latter case, mutual fric- 
tion between the vortex lines and the Ekman layer produce the secondary 
flow, which transports the vortex lines inwards, in agreement with Ref.pQ. 
The analysis agrees qualitatively with the experiments of Tsakadze, but de- 
tailed quantitative agreement remains elusive because the theory does not 
treat the feedback torque on the container self-consistently in the correct 
geometry, nor is it valid in the non-linear regime. 

In this paper, we calculate the spin-down recovery of a superfluid-filled 
container self-consistently and compare our findings with the Tsakadze exper- 
iments [23] • The calculation, whose details can be found in Ref. |281l29j . solves 
the Hall-Vinen-Bekharevich-Khalatnikov (HVBK) two-fluid equations with 
Hall-Vinen mutual friction in arbitrary geometry. The coupled acceleration 
of the fluid and deceleration of the container are followed self-consistently by 
continually updating the fluid boundary conditions in response to the viscous 
back-reaction torque. A summary of the analytic model is given in SJ2J In fJ3] 
the angular velocity versus time, and the dependence of the spin-down time 
on temperature, are compared with measurements from Ref. |25j . Excellent 
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agreement is obtained after correcting for the non-uniform rotation in the 
initial state. 



2 Spin-down theory of a two-component superfluid and its 
container 

Consider a cylindrical or spherical container rotating with angular velocity 
il s about the z axis (unit vector k) . We assume that the superfluid corotates 
uniformly with the container initially (t < 0), as in all such analyses; it turns 
out that experimental data force us to revisit this assumption, as discussed 
in iJ3] At t = 0, the container is spun up impulsively to Q s + SS7 and then left 
to respond to the viscous and frictional torques exerted by the fluid and the 
apparatus. The spin up is linear, with Rossby number e = 5f2/ f2 s <C 1. The 
angular velocity of the container for t > is denoted Q{t) — Q s + 6f2f(t), 
where f(t) is the dimensionless angular velocity perturbation of the con- 
tainer. At time t — 0, when the container is released after being accelerated 
impulsively, the initial velocities of the viscous and inviscid components are 
zero in the rotating frame, and /(0) equals unity. 

In the model, it is assumed that the container initially rotates with Ekman 
number £<1. The Ekman number is defined as 

where v n is the kinematic viscosity the viscous component, and L is the 
appropriate length scale (the radius a for a sphere or the half-height h for a 
cylinder). When E is small, the interior fluid spins up via Ekman pumping. 
It is customary to use the dimensionless variable r, which is the usual time- 
coordinate nor malized by the Ekman time E^ 1 / 2 ^ 1 [101128] . 



2.1 Governing equations 

The response of the superfluid inside the container is described by the incom- 
pressible HVBK equations, which model the superfluid in terms of viscous 
and inviscid components coupled by a mutual friction force. Linearized and 
in the rotating frame, the HVBK equations take the dimensionless form jTTl 

Eg 

E l/2^n + ^ x ^ = + + ^ (2) 



Ot 

d 



v. 



E V'—± + 2k x v s = -Vp s - Pn F , (3) 

OT 

V • v„ = , (4) 

V • v s = . (5) 

where v n)S and p n . s are the bulk velocities and pressures of the viscous (n) 
and inviscid (s) components. We assume laminar flow, where the linearized 
mutual friction force has the Hall-Vinen form [13] 

F = Bkx[kx (v„ - v.)] + B'k x (v„ - v.) . (6) 
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and B, B' are the mutual friction coefficients. The velocity and pressure 
scales are chosen to be L5fi s and pfI s L 2 5fl s respectively, where p is the total 
mass density for the superfluid. The mass densities of the viscous and inviscid 
components, p n and p s , are scaled so that 



To keep the problem analytically tractable, the walls of the container are 
assumed to be smooth, and vortex tension is neglected, an assumption which 
is justified in detail below. The viscous component obeys the usual no-slip 
and no-penetration boundary conditions with respect to the tangential and 
normal velocity respectively. For the inviscid component, no penetration is 
sufficient. 



2.2 Ekman pumping 

In a two-component superfluid, Ekman pumping proceeds as follows. In the 
viscous component, the Coriolis force drives a secondary flow in two layers of 
thickness E 1 ^ 2 L along the upper and lower boundaries of the container. To 
satisfy continuity requirements, the outflow in the boundary layer draws in 
fluid from the interior, spins it up, and cycles it back into the interior on the 
Ekman time-scale E -1 / 2 flj 1 . The inviscid component is spun up by the vis- 
cous component via the mutual friction force. In helium II the mutual friction 
force is strong [B and B' are O(l)], and the inviscid component is dragged 
along with the viscous component as it undergoes Ekman pumping. The az- 
imuthal velocities of both fluid components are "locked together" over the 
Ekman time-scale, behaving essentially as a single fluid [35]. The dimension- 
less quantity that controls the strength of the coupling is 2BE~ 1 / 2 / (2 — B') 
(see below). 

The general analytic solution to has been derived recently 28 , 

generalizing the boundary layer analysis of Greenspan jlOj . The solution is 
most neatly expressed in terms of the total mass flux, defined as 



Pn 



+ Ps 



= I . 



(7) 



v = /0„v n + p s v s . 



(8) 



The azimuthal component of the total mass flux is given by 



v${r, t) 




(9) 
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where we work in cylindrical coordinates denoted (r, 
paper. In ©, we make the auxiliary definitions 



, z) throughout this 



J(r) 



1 

~2 

2BE- 1 ' 2 
2-B' 



I(r) 



h(r) 



I(r) 



h(r) 



pJ{r) 
h(r) 



1/2 



1 



A_(r) [A 2 _(r) + 
' '1 -H 4 (r) 



l 2 



J(r) 



A±(r) 



1 + H 4 (r) 
PnH 2 {r) 



1/2 



[A|(r)-A*(r)] a +4A*(r)A?.(r) 



2pA_(r) [l + iJ 4 (r)] 
x{[A 2 _(r) + A 2 + (r)] [l 



p 1 / 2 
ff(r) 



(B' - 2) 2 



ff 4 (r)] 
- B 2 



(p n B> -2 P y + (p n B) 
2p s B[l + H A {r)] 



-4\ 2 _{r)H 4 (r)} 

1/2 



1/2 



(10) 
(11) 

(12) 
(13) 



H 2 (r) ( Pn B< - 2P) 2 + (p n B)' 



(14) 



Equations ©-(O are valid for a container of arbitrary shape. The fluid 
is contained within the region z — ±h(r). For a sphere, we have h{r) = 
H(r)- 2 = (1 - r 2 ) 1 / 2 . For a cylinder, we have h(r) = H{r) = 1 29!. For 
axisymmetric containers such as spheres and cylinders, the solution is inde- 
pendent of the azimuthal angle </>. 

The solution (JSJ is columnar (i.e., independant of z) as a result of the 
Taylor-Proudman theorem. In Ref. [28 it is demonstrated that for a mutual 
friction force of the form ([6]), each of the fluid components is columnar, 
extending the Taylor-Proudman theorem to the two-fluid problem. 



2.3 Back-reaction torque on container 

Equations (P|)- (IT4l give the solution to ©-© in terms of an arbitrary /(r). 
The motion of the container /(t) is determined by the viscous back reaction 
torque applied by the superfluid. Integrating the viscous stress tensor on the 
boundary z = ±/i(r), torque balance yields 

df(r) 4nK f R a , n ^M, T n 
— — = — / drr h{r) h T , (15) 

where we make the definition 

K= T -f, (16) 
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and R denotes the dimensionless, equatorial radius of the vessel (a/h for a 
cylinder of height 2h and radius a; unity for a sphere). The constant fric- 
tional torque in the apparatus is denoted by To in dimensionless units and is 
negative |29] . // and I c denote the dimensionless nominal moments of inertia 
of the fluid (rotating as if it were a rigid body of uniform density p) and the 
container respectively, scaled to pL 5 . For a sphere, we have // = 87r/15; for 
a cylinder, we have If — ttR. 

Equations ([9]) and (f!5|) constitute a closed pair of integro-differential equa- 
tions for the unknowns v$(r,r) and /(t). Substituting (fT4"]) into ([15]). the 
integral equation of motion for the container is [H] 

/ (r) = -K f dr' [g A (r - r') + g B (r - r')] f(r') + T r + 1 , (17) 
Jo 

with 

4tt drr 3 J(r) [e w +W r -e^-W 7 "] 

9 ( T ) = ~t~ / n n ' ( 18 ) 

If Jo w+(r) -w_(r) 

B 4tt /•« drr 3 {^_(r) [e^K _ i] _ ^ +(r ) [ e ^(» - l] } 



If Jo w+(r) -w_(r) 

An overdot symbolizes a derivative with respect to r. 

Equation (|17l) captures self-consistently the back-reaction of the super- 
fluid on the container and vice versa, e.g. by continuously varying the bound- 
ary conditions on the superfiuid as fi(r) evolves. It is straightforward to solve 
(fT7)l numerically, by guessing an initial trial function for /(r), substituting it 
into the right-hand side of (IT71) . updating /(r) via underrelaxation [TB], and 
iterating. This procedure is described in Ref. [25] . 



2.4 Vortex tension 

In deriving (|17[) . the vortex tension has been neglected. The reasons for this 
are discussed in detail in Ref. [28], and we re- iterate them here. 

In the solution of Reisenegger [Tj[5] for slow Ekman pumping between 
slowly accelerating plates, the tension coefficient v s only appears in the in- 
viscid velocity component in the boundary layer 0. The qualitative features 
of the slow-acceleration boundary problem are expected to carry over to the 
impulsive spin- up problem, just like they do in the classical case [3]. Hence, 
v 8 does not appear in the solution for the interior flow or the viscous velocity 
component in the boundary layer. Nor, consequently, does it appear in the 
spin- up time. [In section 2 of Ref. [55], the results of Reisenegger [B5] are 

1 In fact, v s only enters in the combination v s /v n , which is the ratio of the tension 
force to the viscous force. Although the length-scale of the system may enter the 
force ratio in principle, it does not do so in practice. Chandler and Baym ^ wrote 
the linearized tension force as 2Qv s {d / 'dz) 2 £ (where £ is the vortex displacement 
vector) , which is proportional to inverse length squared and hence scales identically 
with the viscous force. Alternatively, the Hall-Vinen form of the tension is given 
by v s ujs x (Vxii), which also scales in inverse proportion to length squared. 
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recovered despite assuming v s = 0.] Thus, vortex tension is not neglected on 
the basis of scaling but by realizing that, during the spin up of smooth-walled 
containers, the tension v s appears only in the inviscid component boundary 
layer and therefore does not affect the gross dynamics (e.g. spin-up time) q 

Also, because the torque is calculated from the viscous-component bound- 
ary layer (which is independent of u s ), the feedback calculated in 32.31 and 
Ref . [21] is also independent of v s . In §S\ we find excellent agreement between 
theory and experiment, supporting the hypothesis that vortex tension does 
not play a significant role in the linear spin-up of superfluids in smooth- walled 
containers at temperatures below 1.6 K. 

Nevertheless, it is interesting to evaluate the characteristic strength of the 
tension force relative to other forces. For a neutron superfluid we have v s ~ 
lCP'm 2 s _1 , therefore, v s jv n ~ 10~ 12 . Clearly, for astrophysical systems like 
neutron stars, we find that the vortex tension is negligible. However, for 
laboratory experiments with helium II, we find 0.5 < v s jv n < 6.1 for 1.3 K < 
T < 2.1 K. Hence, if it were not for the form of the tension force discussed in 
the previous paragraph, it would have to be included. An interesting avenue 
for future work is to repeat the calculation in Ref. [35] keeping the vortex 
tension, to clarify its role in impulsive spin up. 

Because vortex lines must emerge perpendicular to the surface, there 
is a region close to the surface where vortex lines are not vertical in non- 
parallel-plate geometries. The thickness of this region is governed by the 
vortex line tension. By analogy with the classical viscous penetration depth 
fin = {v n / 'ft) 1 / 2 , we define a vortex tension penetration depth 5 S , which in 
helium II takes the value 

in agreement with Ref. [115]. Because v s and v n are comparable in helium II, 
this region is of approximately the same thickness as the Ekman layer. 



3 Helium II experiments 

3.1 Vessel rotation curve 

In the Tsakadze experiments, spherical and cylindrical containers are studied 
under a variety of initial conditions. Only those experiments in which the 
container is initially rotating are examined in this paper, to ensure that E 
and e are small and the analytic solution in ^\ is valid. For the helium II 
experiments of the Tsakadze group [25) . we find 

£ = 3.5x 10- 5 ( ) ( ° \ 1 [±- V 2 , (21) 

V5.5 x 10- 8 Pasy V lrads / \ 4cn V 

2 The ratio of the tension force to the mutual friction force is also small (of order 
v a /L 2 Q ~ E for the linear problem in helium II). 
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at T = 1.6 K. For 1.2K <T< 1.6 K, one measures 1.197 x 1CT 8 < v n < 
8.997 x 10~ 9 using vibrating wire viscometers [22,[5] (see also Donnelly's 
websitcQ). Therefore the requirement E <C 1 is valid for the experiments we 
consider here. 




100 200 300 400 100 200 300 400 



t (s) t (s) 

Fig. 1 Angular velocity fl(t) (in rads^ 1 ) of a spherical container versus time (in s) 
following an impulsive acceleration from fi a — 3.0rads _1 to Q s + 8Q = 4.2rads _1 
at t = 0. The dots and crosses are the experimental data from Ref. [25] for an empty 
vessel and one filled with helium II at T — 1.57 K. In both cases, experimental error 
bars are smaller than the plotted symbols. The heavy solid curve is the theoretical 
solution presented in Sj2] (Left panel.) Superfluid coefficients taken from Table 1. 
(Right panel.) As for left panel, but with B = 1.22 x 1CT 2 and B' = 2.31 x 10" 3 . 



In the left panel of Figure 1, we plot the spin-down recovery of a spher- 
ical container of radius a = 1.7 cm, filled with helium II at 1.57 K, after an 
impulsive increase in angular velocity from Q s = 3.0rads _1 to fl s + 5Q = 
4.2rads _1 . The dots and crosses represent experimental data for the empty 
vessel and helium-filled vessel respectively, corresponding to curves 3 and 2 in 
Figure 7 of Ref. [35] . The slope of the spin-down curve for the empty container 
as a result of friction in the apparatus is &QQ s E x I 2 Tq = —0.002 rad s~ 2 . 
The theoretical spin down for a spherical container calculated from (|17p is 
over plotted as a heavy solid curve, using the values of the superfluid coef- 
ficients quoted in Table 1. The crosses and solid curve approach a constant 
slope fl = 5f2f2 s E 1 / 2 To(l + K)~ x asymptotically. Comparing the asymptotic 
slopes of the crosses and dots gives K = 0.8. 



Table 1 Superfluid coefficients for helium II at T = 1.57 K 

Quantity Value Units Source 

p„ 21.2 kgm~ 3 Mavnard et al. [14] 

p 145.2 kgm" 3 Mavnard et al. [14] 

B 1.222 Barenghi & Donnelly [2] 

B' 0.231 Barenghi & Donnelly [2] 

■q 1.26 x 10" 6 Pas Tough et al. [21] 



3 http://darkwing.uoregon.edu/ rjd/vaporl3.htm 
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The left panel of Figure 1 demonstrates qualitative agreement between 
theory and experiment. The heavy solid curve parallels the crosses asymp- 
totically, but the vertical intercepts extrapolated back from the asymptotic 
solutions disagree by 0.32 rads -2 . It is temping to speculate that the dis- 
crepancy arises because we are adopting inappropriate values of the hard-to- 
measure transport coefficients r), B and B' . For example, B and B' decrease 
by several orders of magnitude when the Donnelly- Glaberson instability [7] 
is excited by a counterflow and a vortex tangle forms p2B] . 

To test this idea, we reduce B and B' 10 -2 -fold and re-plot the data 
and theory. The agreement with experiment worsens: the heavy curve now 
takes longer to reach its constant- J? asymptote, and the vertical intercept of 
the back-extrapolated asymptote is 0.1 rads -1 lower than in the left panel. 
Exhaustive trials reveal that the experimental data cannot be matched by 
changing the superfluid coefficients. 

Another explanation for the discrepancy is that the initial conditions are 
incorrect. An important clue that this is in fact the case follows from close 
inspection of the experimental data. When the angular impulse removed by 
the frictional torque is subtracted, the final angular momentum of the fluid 
and container is much less than the initial angular momentum, given fl s = 
3.0rads -1 . It is reported in Ref. [25] that J? rises to its maximum at t ~ 10s, 
whereas the Ekman time is ~ 40 s. Hence the fluid partially spins up before 
the container is released, partially invalidating the impulsive initial conditions 
assumed in the theory in Because the frictional torque is known, the 
corrected initial angular velocity of the superfluid can be calculated from 
angular momentum conservation 0. 




t(s) 



Fig. 2 Angular velocity of a spherical container filled with helium II (at T — 
1.57 K) as a function of time, with initial conditions (O s — 3.68 rads -1 and fi s + 
SO = 4.25 rads -1 ) adjusted to account for the incomplete pre-release spin up. The 
blac k cross marks the point on the curve where the spin-up time quoted in Ref. 
[25] is read off. The theoretical solution asymptotes to the dotted line. 



4 The 'initial' (at r = when the container is released after it is accelerated) 
angular momentum of the vessel (known) plus fluid (unknown) equals the 'final' 
(at any instant tj in the asymptotic regime) angular momentum of the vessel 
[kno wn from Q(t;)] plus fluid [known in terms of fi(ft) from the theory in [J2]and 
[29| in the asymptotic regime] plus the angular impulse taken away by the friction 
(Tot/). As the angular impulse is given self-consistently and uniquely by the theory, 
subtracting it to find the initial angular momentum of the fluid is equivalent to 
adjusting the single unknown O a until theory and data agree. 
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In Figure 2, we re-plot the solid black curve in Figure 1 with modified 
initial conditions for the superfluid. We find that values of Q s = 3.68 rads -1 
and fl s + 8fl = 4.25 rads -1 give near perfect agreement (error w 0.5%) be- 
tween theory and experiment. In this situation, the superfluid and container 
rotate at 3.0 rads -1 before spinning up. At time r = 0, when the container 
is released after being accelerated, its angular velocity is 4.25 rads -1 , while 
the angular velocity of the superfluid has increased to 3.68 rads -1 . The close 
agreement obtained between theory and experiment in the left panel of Fig- 
ure 2 is strong evidence in support of the hypothesis that helium II spins up 
via laminar Ekman pumping. When the viscous component density fraction 
is small (e.g. at T — 1.57 K), the spin- up flow is dominated by Ekman- like 
circulation in the inviscid component [H] . The result also suggests that vor- 
tex pinning and tension (which are neglected in our theoretical analysis) are 
not needed to reproduce the experimental data in Ref. |25] . 




t(s) 



Fig. 3 Same as Figure 2, but for a cylinder with h — a — 1.7 cm and S7 S = 
3.8 rads -1 . The agreement between theory and data is poorer than for a sphere. 

We note in passing that the shape of the vessel must be treated correctly, if 
the theory is to reproduce the data. Figure 3 displays the spin down of a cylin- 
der filled with helium II at 1.57 K, which is identical to the apparatus modeled 
in Figure 2 (with half- height and cylindrical radius h = a = 1.7 cm) except 
for its shape. Even after correcting the initial conditions (f2 s = 3.8 rads -1 
and S7 s +5f2 = 4.25 rads -1 ), the theory agrees poorly (maximum error « 3%) 
with the experimental data. Spheres and cylinders produce characteristically 
different spin-down curves, because the torque varies with latitude in a sphere 

mi- 



3.2 Superfluid spin-up time 

The measured spin-up times to for experiments at six different temperatures 
(with all other experimental parameters identical to Figure 2) are plotted as 
points with error bars in Figure 4. The data are taken from curve 2 of Figure 
9 in Ref. [55] • The definition of to adopted in Ref. [25] refers to the point at 
which the spin-down curve reaches its asymptotic regime, drawn as a dotted 
line in Figure 2. Clearly, there is scope for ambiguity in identifying exactly 
when this occurs. In Ref. [25] and Figure 4, the value of to — 119.65 s quoted 
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at T = 1.57 K corresponds to the point marked with a cross in Figure 2, 
where the difference between J?(i) and the long-term linear spin down drops 
to 11% of its original value. To remain consistent, we calibrate against this 
point and adopt the foregoing definition of to throughout the rest of this 
paper. 




1.4 1.6 1.8 2.0 2.2 



T(K) 

Fig. 4 Spin- up time to (in s) versus temperature, T (in K). Experimental mea- 
surements from Ref. 1251 are plotted as black points with error bars. The spin-up 
time calculated from (|17p and calibrated against the cross in Figure 2 is plotted as 
a heavy black curve. 

In Figure 4, the theoretical spin-up time [calculated from (fTT|) ] is plot- 
ted versus temperature as a heavy black curve. Good agreement with the 
data is obtained for temperatures between 1.46 K < T < 1.8 K, where the 
inviscid component density fraction is less than 0.3. Above 2.0 K, the viscous 
component dominates and the agreement worsens. Turbulent drag may be 
responsible. In a classical fluid, the Ekman layer becomes unstable when the 
Ekman number drops below the critical value [5111 9j 

E = e 2 (56.3 + 58.4er 2 . (22) 

When the flow becomes turbulent, the drag is reduced, lengthening to. Above 
the critical temperature, where the flow is classical, the condition (|2"2"|) applies 
directly. Below the critical temperature, transition to turbulence in superflu- 
ids is not well understood, although it typically occurs at Ekman numbers 
similar to (or even higher than) the classical critical value [6lf2Tj . Reisenegger 
[T5] suggested a critical value of E w 2.5 x 10~ 3 e 2 , reasoning that the vortex 
lines should be approximately parallel to the rotation axis in the laminar 
regime, although it should be noted that this is never strictly allowed by the 
boundary conditions in a sphere (where the vortex lines emerge perpendic- 
ular to the surface). Under the experimental conditions in Figures 2 and 4, 
the above results imply that the flow should be turbulent for T > 1.6 K. 

For the spin down of smooth-walled spherical containers, Tsakadze and 
Tsakadze |25j derived the empirical relation 

n >*= A (^^)\j) ° h (l + ^), (23) 

where m is the mass of a helium atom, and h is Planck's constant divided 
by 27T. The remaining constants are fitted by the method of least squares, 



12 



yielding A = 1.0±0.1, 7 = 0.40±0.05, a = 0.25±0.01 and c = SiO^rad" 1 s. 
In the analytic theory presented in inland Ref. [251129] . to scales as the Ekman 
time, implying ft s to ~ (a 2 /2s) 1 / 2 , in close agreement with (1231) . For linear 
spin- up, to does not depend on 6 ft. However, in deriving (f2"3")l . flows with e 
as high as 0.3 were studied, where the linearity assumption is questionable, 
the flow may be turbulent, and to does depend on 6 ft. 

The temperature dependence in (l2"3")l is contained in p„/p. In the left 




1.0 1.2 1.4 1.6 1.8 2.0 1.0 1.2 1.4 1.6 1.8 2.0 



(P„/P)- - 25 (P„/P)-"- 25 

Fig. 5 Spin-up time to (in s) versus the reciprocal of the fourth root of the density 
fraction of the viscous component, (p n /p)~ ' ■ Left panel: fl s — 3.0rads _1 , Q s + 
8fl = 4.2rads~ 1 (as in Figure 3). Right panel: fl s = 5.0rads _1 , fl s + SO = 
6.5rads _1 . The dashed curves represent t he e mpirical relation (|23p . Solid curves 
are the theoretical results calculated from (|17|l . 



panel of Figure 5, values from the horizontal axis in Figure 4 are converted 
from T to (p n / p)~ - 25 , using generally accepted values of p(T) and p n {T), 
e.g. contained in Ref. [14]. The experimental data are plotted as points with 
error bars, the theoretical curve [calculated from ([T7|] is plotted as a heavy 
black curve, and the empirical relation (1231) is over plotted as a dashed curve. 
Note that temperature now increases from right to left. The critical value 
£ k 2.5 x l(F 3 e 2 occurs at (/0,i/p>)~ ' 25 = 1.58, below which turbulence is 
expected. Equation (1T7|) describes the data better than equation (12"3"|) for 
these parameters, suggesting that (|2"31 is unreliable when extrapolating to 
parameters outside those for which it is derived. 

In the right panel of Figure 5, we compare the theoretical predictions from 
([TTf with the specific experiments used to fit A, 7, a and c in (|2"3"]l . A sphere of 
radius a — 1.7 cm is spun up from ft s = 5.0rads _1 to ft s + 6 ft = 6.5rads _1 . 
The experimental data from Figure 10 in Ref. [25 are plotted as points 
with error bars. In this case, the theory departs noticeably from the data. 
However, for these parameters the threshold E « 2.5 x 10 _3 e 2 is achieved 
at T = 1.25 K and hence {p n /p)~ - 25 = 2.32. This suggests that all the 
data points in the right panel of Figure 4, and therefore A, 7, a and c in the 
empirical relation (j23[) . apply to the turbulent regime. Unfortunately, because 
there are no measurements in the laminar regime [(p n / p) ' 25 > 2.32], the 
theory cannot be conclusively tested with these data. 

We close by remarking on the radius of the vessel used in the experiments. 
In Ref. [25], a value of 3.4 cm is quoted as both the radius and the diameter 
[19] . In this paper we assume a — 1.7cm everywhere. The larger value of 
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a = 3.4 cm leads to double the spin-up time, making agreement between 
theory and experiment very difficult to obtain. Under this interpretation, 
the radius quoted in equation (6) of Ref. [25] refers instead to the diameter. 
An appropriate correction is made in equation (f2"3")l of this paper. 



4 Conclusions 

A self-consistent analytic solution of the HVBK equations for the spin up 
of helium II enclosed in a vessel of arbitrary geometry, including the vis- 
cous feedback torque, is compared with experimental data from the classic 
experiments conducted by Tsakadze and Tsakadze [25] . Excellent agreement 
between the theoretical and observed spin-down of the container is obtained 
at T = 1.57 K. The agreement supports the hypothesis that laminar Ekman 
pumping occurs, which vortex tension is negligible and the container walls 
are smooth enough to avoid much pinning. Because the container initially 
spins up over ~ 10 s, a fair fraction of the Ekman time (~ 40 s), the acceler- 
ation event is not truly impulsive, and the theoretical initial conditions must 
be corrected by angular momentum conservation arguments |29] . 

As one approaches the superfluid transition temperature, the theoretical 
and observed spin-down times diverge. It is suggested that this discrepancy 
is caused by the onset of turbulence as r\ and p n /p increase, consistent with 
a critical Ekman number E rj 2.5 x 10~ 3 e 2 proposed in Ref. [T!5]. Several 
observable differences are present if the flow is turbulent: (1) the spin-up 
time depends on Sf2, and (2) the sphere and cylinder have similar spin- 
down curves, because latitudinal variation of the hydrodynamic torque in 
Ekman pumping is wiped out to some extent. Further experiments, conducted 
specifically in the range e < E 1 ! 2 <C 1 are required to test the theory fully. 
The coupled spin up of a superfluid and its container is a vital component 
in understanding the nature of vortex pinning and nucleation in terrestrial 
settings and the recovery process following rotational glitches in neutron 
stars. 
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